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Abstract: Characteristics of uncertainty and imprecision, even imperfection is 
presented from knowledge acquisition in map reconstruction of autonomous mobile 
robots. Especially in the course of building grid map using sonar, this characteristic 
of uncertainty is especially servere. Jean Dezert and Florentin Smarandache have 
recently proposed a new information fusion theory (DSmT), whose greatest merit is 
to deal with uncertainty and conflict of information, and also proposed a series of 
proportional conflict redistribution rules (PRCI~PCRS5), therein, presently PCR5 
is the most precise rule to deal with conflict factor according to its authors, though 
the complexity of computation might be increased correspondingly. In this chapter, 
according to the fusion machine based on the theory of DSmT coupled with PCR5, 
we not only can fuse information of the same reliable degree from homogeneous 
or heterogeneous sensors, but also the different reliable degree of evidential sources 
with the discounting theory. Then we established the belief model for sonar grid map, 
and constructed the generalized basic belief assignment function (gbbaf). Pioneer IT 
virtual mobile robot with 16 sonar range finders on itself served as the experiment 
platform, which evolves in a virtual environment with some obstacles (discernable 
objects) and 8D Map was rebuilt online with our self-developing software platform. 
At the same time, we also compare it from other methods (i.e. Probability theory, 
Fuzzy theory and Dempster-Shafer Theory (DST)). The results of the comparison 
shows the new tool to have a better performance in map reconstruction of mobile 
robot. It also supplied with a foundation to study the Self-Localization And Mapping 
(SLAM) problem with the new tool further. 
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14.1 Introduction 


The study on exploration of entirely unknown environment for intelligent mobile robots has been 
a popular and difficult subject for experts in the robotic field for a long time. Robots do not 
know the environment around themselves, that is, they have no experienced knowledge about the 
environment such as size, shape, layout of the environment, and also no signs such as beacons, 
landmarks, allowing them to determine their location about robot within the environment. 
Thus, the relation between self-localization and map building for mobile robot is like the chicken 
and egg problem [3, 16]. This is because if the mobile robot builds the map of the environment, 
then it must know the real position of itself within the environment; at the same time, if the 
robot wants to know its own position, then it must have a referenced map of the environment. 
Though it is hard to answer this question, some intelligent sensors such as odometer, electronic 
compass, sonar detector, laser range finder and vision sensor are installed on the mobile robot 
as if a person has perceptive organs. 

How to manage and utilize this perceptive information acquired by organs, it’s a new subject 
in information fusion, which will play an important role herein. As far as we know, experts 
have not yet given a unified expression. Just aiming to the practical field or system, proposed 
architecture of control such as hierarchical, concentrative, distributive and composite, and then 
according to the different integrated hierarchy, we compared the validity of all kinds of classical 
(Probability) and intelligent (Fuzzy, Neural-Networks (NN), Rough Set theory, Dempster-Shafer 
theory (DST), etc.) arithmetic. As far as the mobile robot is concerned, the popular arithmetic 
of self-localization in an unknown environment relying on interoceptive sensors (odometer, elec- 
tronic compass) and exteroceptive sensors (sonar detector, laser range finder and visual sensor) 
is Markov location [10] or Monte Carlo location [28]. The map of the environment is built by 
applying some arithmetic such as Probability theory, Fuzzy Set theory and DST. The informa- 
tion of environment can be expressed as grid map, geometrical feature or topological map, etc., 
where the grid map is the most popular arithmetic expression [8, 9]. In this chapter, a new tool 
of the Fusion Machine based on DSmT [5, 6, 22] coupling with PCR5 is introduced to apply 
to the map reconstruction of mobile robots. DSmT mentioned here that has been proposed 
by Jean Dezert and Florentin Smarandache based on Bayesian theory and Dempster-Shafer 
theory [21] recently is a general, flexible and valid arithmetic of fusion. Its largest advantage 
is that it can deal with uncertain and imprecise information effectively, which supplies with a 
powerful tool to deal with uncertain information acquired by sonar detector in the course of 
building the grid map. Moreover, through the rule of PCR5, which is also proposed by Jean 
Dezert and Florentin Smarandache [23-25], we can refine and redistribute the conflict mass to 
improve the precision and correctness of fusion. The comparison of the new tool from other 
methods is done to testify it to have a better performance to solve the puzzle. 


14.2 The fusion machine 


14.2.1 General principle 


At first, here the fusion machine is referred to a theory tool to combine and integrate the 
imperfect information without preprocessing it (ie. filter the information) according to the 
different combination rules (i.e. DST, DSmT, etc.). It even redistributes the conflict masses to 
other basic belief masses according to the constraints of system using the different redistribution 
rules (i.e. PCRI~PCRS5, minC [2], WAO [11], etc.). Of course, how to adopt the fusion rule must 
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be considered according to the different application. Here we consider the application, and give a 
special fusion machine (shown in Fig. 14.1). In Fig. 14.1, k sources of evidences (i.e. the inputs) 
provide basic belief assignments over a propositional space generated by elements of a frame 
of discernment and set operators endowed with eventually a given set of integrity constraints, 
which depend on the nature of elements of the frame. The set of belief assignments need then to 
be combined with a fusion operator. Since in general the combination of uncertain information 
yields a degree of conflict, say K, between sources, this conflict must be managed by the fusion 
operator/machine. The way the conflict is managed is the key of the fusion step and makes 
the difference between the fusion machines. The fusion can be performed globally /optimally 
(when combining the sources in one derivation step all together) or sequentially (one source after 
another as in Fig. 14.1). The sequential fusion processing (well adapted for temporal fusion) is 
natural and more simple than the global fusion but in general remains only suboptimal if the 
fusion rule chosen is not associative, which is the case for most of fusion rules, but Dempster’s 
rule. In this chapter, the sequential fusion based on the PCR5 rule is chosen because PCR5 has 
shown good performances in works and because the sequential fusion is much more simple to 
implement and to test. The optimal (global) PCR5 fusion rule formula for k sources is possible 
and has also been proposed [23] but is much more difficult to implement and has not been 
tested yet. A more efficient PCR rule (denoted PCR6) proposed very recently by Martin and 
Osswald in [15], which outperforms PCR5, could be advantageously used in the fusion machine 
instead PCR5. Such idea is currently under investigation and new results will be reported in 
a forthcoming publication. We present in more details in next section the DSmT-based fusion 
machine. 


Figure 14.1: A kind of sequential fusion machine 


14.2.2 Basis of DSmT 


DSmT (Dezert-Smarandache Theory) is a new, general and flexible arithmetic of fusion, which 
can solve the fusion problem of different tiers including data-tier, feature-tier and decision-tier, 
and even, not only can solve the static problem of fusion, but also can solve the dynamic one. 
Especially, it has a prominent merit that it can deal with uncertain and highly conflicting 
information [5, 6, 22]. 
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14.2.2.1 Simple review of DSmT 


1) Let O = {0),02,--- ,A,}, be the frame of discernment, which includes n finite focal elements 
6;(i = 1,--- ,n). Because the focal elements are not precisely defined and separated, so that no 
refinement of © in a new larger set (,¢f of disjoint elementary hypotheses is possible. 

2) The hyper-power set D® is defined as the set of all compositions built from elements of 
with U and (© generates D®© under operators U and M) operators such that 


a) 0,01,02,-°° On E D®. 
b) If A,B € D®, then ANB € D® and AUBE D®. 
c) No other elements belong to D®, except those obtained by using rules a) or b). 


3) General belief and plausibility functions 

Let O = {0),42,--- ,An} be the general frame of discernment. For every evidential source 
S, let us define a set of map of m(-) : D° — [0,1] associated to it (abandoning Shafer’s model) 
by assuming here that the fuzzy/vague/relative nature of elements 6;(i = 1,--- ,) can be non- 
exclusive, as well as no refinement of © into a new finer exclusive frame of discernment ©, is 
possible. The mapping m(-) is called a generalized basic belief assignment function if it satisfies 


m(0)=0 = and iS m(A) = 1, 


AED® 


then m(A) is called A’s generalized basic belief assignment function (gbbaf). The general belief 
function and plausibility function are defined respectively in almost the same manner as within 
the DST, ice. 

Bel(A)= S>  m(B) (14.1) 


BeED®,BCA 


PY(A)= S> mB) (14.2) 


BeED®,BNAFAO 


4) Classical (free) DSm rule of combination 
Let M/(®) be a free DSm model. The classical (free) DSm rule of combination (denoted 
(DSmC) for short) for k > 2 sources is given VA 4 0, andA € D®as follows: 


k 
mys(o)(A) = fm @---@mel(A)= SS [mil Xi) (14.3) 
X15 ,XpED2 1=1 
X{N AX pHA 


14.2.2.2 Fusion of unreliable sources 


1) On the necessity of discounting sources 

In fact, sources of information are unreliable in real systems due to the sources with different 
knowledge and experience. For example, from the point of view of the mobile robots’ sensors, 
the metrical precision and resolution with laser range finder are both higher than that with 
sonar sensor. Even if they are the same sonar sensors, then they have also different precision 
due to the manufacturing and other factors. Under this condition, if we treat data of unreliable 
information sources as data of reliable sources to be fused, then the result is very unreliable 


14.2. THE FUSION MACHINE 347 


and even reverses. Thus, unreliable resources must be considered, and then DSmT based on 
the discounting method [7, 12, 21, 26] does well in dealing with unreliable sensors. 
2) Principle of discounting method 

Let’s consider k evidential sources of information ($1, S2,...,S%), here we work out a uniform 
way in dealing with the homogeneous and heterogeneous information sources. So we get the 
discernment frame © = {61,62,--» ,On}, m/(-) is the basic belief assignment, let m;(-) : D9 — 
[0,1] be a set of maps, and let p; represent reliable degree supported by S; (¢ = 1,2,...,k), 
considering }7 yepe mi(A) = 1, let I, = 0; U02U---U 6, express the total ignorance, and then 
let m?(I:) = 1 — p; + pimi(Iz) represent the belief assignment of the total ignorance for global 
system (after discounting), and then this is because of existing occurrence of malfunction, that 
is, )) 4epe Mi(A) = pi, we assign the quantity 1—>p; to the total ignorance again. Thus, the rule 
of combination for DSmT based on discounting method with k > 2 evidential sources is given as 


in the formula (14.3), i.e. the conjunctive consensus on the hyper-power set by ms #(@) 9) =0 
and VA £0 € D®, 
k 
mis@(A) =lmfe...emi(A4)= So [piri Xi) (14.4) 


Xyy°,.XpED2 i=1 
Xp NX p=A 


14.2.3. The PCR5 fusion rule 


When integrity constraints are introduced in the model, one has to deal with the conflicting 
masses, i.e. all the masses that would become assigned to the empty set through the DSmC 
rule. Many fusion rules (mostly based on Shafer’s model) have been proposed [20] for managing 
the conflict. Among these rules, Dempster’s rule [21] redistributes the total conflicting mass 
over all propositions of 2° through a simple normalization step. This rule has been the source 
of debates and criticisms because of its unexpected/counter-intuitive behavior in some cases. 
Many alternatives have then been proposed [20, 22] for overcoming this drawback. In DSmT, 
we have first extended the Dubois & Prade’s rule [7, 22] for taking into account any integrity 
constraints in the model and also the possible dynamicity of the model and the frame. This first 
general fusion rule, called DSmH (DSm Hybrid) rule, consists just in transferring the partial 
conflicts onto the partial ignorances!. The DSmH rule has been recently and advantageously 
replaced by the more sophisticated Proportional Conflict Redistribution rule no.5 (PCRS). 
According to Smarandache and Dezert, PCR5 does a better redistribution of the conflicting 
mass than Dempster’s rule since PCR5 goes backwards on the tracks of the conjunctive rule 
and redistributes the partial conflicting masses only to the sets involved in the conflict and 
proportionally to their masses put in the conflict, considering the conjunctive normal form of 
the partial conflict. PCR5 is quasi-associative and preserves the neutral impact of the vacuous 
belief assignment. Since PCRS is presented in details in [23], we just remind PCR5 rule for only 
two sources”: mpcrs(0) =0 , and for all X € G \ {0}, 


m3(X)mo(¥)_, _m3(X)ri(¥) 


mpcrs(X) = m12(X) + m1(X) + m2(Y) m2(X) +my(Y) 


(14.5) 


YEG\{x} | 
XNY=0 


Partial ignorance being the disjunction of elements involved in the partial conflicts. 
2A general expression of PCR5 for an arbitrary number (s > 2) of sources can be found in [23]. 
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where all sets are in canonical form and m2(X) = 0 x,,x,ece ™1(X1) - M2(X2) corresponds 
X1AX5 
to the conjunctive consensus on X between the two sources and where all denominators are 


different from zero. If a denominator is zero, that fraction is discarded. 


Figure 14.2: Sketch of the principle of sonar 


14.3. Modeling of Sonar Grid Map Building Based on DSmT 


Here we mainly discuss a sonar sensor, whose working principle (shown as Fig. 14.2) is: pro- 
ducing sheaves of cone-shaped wave and detecting the objects by receiving the reflected wave. 
Due to the restriction of sonar physical characteristic, metrical data has uncertainty as follows: 
a) Beside its own error of making, the influence of external environment is also very great, for 
example, temperature, humidity, atmospheric pressure and so on. 

b) Because the sound wave spreads outwards through a form of loudspeaker, and there exists 
a cone-shaped angle, we cannot know the true position of object detected among the fan-shaped 
area, with the enlargement of distance between sonar and it. 

c) The use of many sonar sensors will result in interference with each other. For example, 
when the i-th sonar gives out detecting wave towards an object of irregular shape, if the angle 
of incidence is too large, the sonar wave might be reflected out of the receiving range of the 7-th 
sonar sensor or also might be received by other sonar sensors. 

d) Because sonar sensors utilize the reflection principle of sound wave, if the object absorbs 
most of heavy sound wave, the sonar sensor might be invalid. 


Pointing to the characteristics of sonar’s measurement, we construct a model of uncertain 
information acquired from grid map using sonar based on DSmT. Here we suppose there are 
two focal elements in system, that is, 0 = {01,02}. Where 6; means grid is empty, 62 means 
occupied, and then we can get its hyper-power set D° = {0,01 9 02,01, 02,01 U 02}. Every grid 
in environment is scanned k > 5 times, each of which is viewed as source of evidence. Then 
we may define a set of map aiming to every source of evidence and construct the general basic 
belief assignment functions (gbbaf) as follows: m(61) is defined as the gbbaf for grid-unoccupied 
(empty); m(@2) is defined as the gbbaf for grid-occupied; m(61M 62) is defined as the gbbaf for 
holding grid-unoccupied and occupied simultaneous (conflict). m(61 U 42) is defined as the gb- 
baf for grid-ignorance due to the restriction of knowledge and present experience (here referring 
to the gbbaf for these grids still not scanned presently), it reflects the degree of ignorance of 
egrid-unoccupied or occupied. 


The gbbaf of a set of map m(-) : D° — [0,1] is constructed by authors such as the formulae 
(14.6)~(14.9) according to sonar physical characteristics. 


14.3. MODELING OF SONAR GRID MAP BUILDING BASED ON DSMT 


Rmin <p <R< Rmax 
1—(p/R)?)d if 
(1=(p/RP)A i cere 


0 otherwise 


eo 3Pv(p—R)? y if Rmin < p < R +e < Rmax 
0<@<wu/2 


0 otherwise 


m<xp<R< 
fi pies : te SpsRkS Raz 


m(0, 7 62) = 0<60<w/2 
0 otherwise 
R< < Rmax 
tanh(2(p — R))A if ae 
m(6, U 02) = 0<0<u/2 


0 otherwise 


where \ = E(0) = O(6) is given by (see [8] for justification) 


ja 1—(20/w)? if 0<|6|<w/2 
~ 10. otherwise 
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(14.6) 


(14.7) 


(14.8) 


(14.9) 


(14.10) 


where py in formula (14.7) is defined as an environment adjusting variable, that is, the less the 
object is in environment, the greater the variable p, is, and makes the function of m(62) more 
sensitive. Here let p, be one. E(-) and O(-) are expressed as the Effect Function of p, 6 to grid’s 
empty or occupancy. In order to insure the sum of all masses to be one, we must renormalize it. 
The analysis on the characteristics of gbbaf are shown as Fig. 14.3~Fig. 14.7, when R = 1.5m. 


Ll L L 
0 0.5 1 1.5 2 25 3 


Figure 14.3: m(61) as function of p given by (14.6) 
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L L L L 
0 0.5 H| 1.5 2 25 3 


Figure 14.5: m(61M @2) as function of p given by (14.8) 


Seen from Fig. 14.3, m(0,) has a falling tendency with the increasing of distance between 
grid and sonar, and has the maximum at R,,j, and zero at R. From the point of view of the 
working principle of sonar, the more the distance between them approaches the measured value, 
the more that grid might be occupied. Thus the probability that grid indicated is empty is very 
low, of course the gbbaf of grid-unoccupied is given a low value. 


From Fig. 14.4, m(@2) takes on the distribution of Gaussian function with respect to the 
addition of distance between them, has the maximum at R, which answers for the characteristic 
of sonar acquiring information. 


From Fig. 14.5, m(@, 9 @2) takes on the distribution of a parabola function with respect 
to the addition of distance between them. In fact, when m(@;) equals m(62), m(@1 9 62) has 
the maximum there. But it is very difficult and unnecessary to find the point of intersection 
of the two functions. Generally, we let the position of R — 2¢ replace the point of intersection. 
Experience indicates that its approximate value is more rational. 
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L L L L 
0 0.5 1 15 2 2.5 3 


Figure 14.6: m(6, U @2) as function of p given by (14.9) 


0 Ll L L i L 
-15 -10 -5 0 5 10 15 


Figure 14.7: as function of @ given by (14.10) 


From Fig. 14.6, m(0,U@2) takes on the distribution of the hyperbola function with respect to 
the addition of distance between them, and zero at R. This function reflects well the ignorance of 
grid information at R < p < Rmaz. The relation between @ and 4) is reflected in Fig. 14.7, where 
the more the position of grid approaches the central axis, the greater \ becomes, that is, the 
greater the contribution to belief assignment is. Otherwise, the lower it is. In short, the general 
basic belief assignment functions (gbbaf) entirely fit with the characteristic of sonar acquiring 
information. This supplies a theoretic foundation for dealing with uncertain information in grid 
map building. 


14.4 Sonar Grid Map Building Based on Other Methods 


To apply the probability theory and fuzzy set theory to map building, at first, two functions of 
uncertainty are introduced. Here the working environment U of robot is separated into m x n 
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Figure 14.8: Sketch of the layout of sonars 


rectangle grids of same size. Every grid is represented by Gi;, U = {Gi;|¢ € [1, mJ], 7 € [1, n]}, 
according to the reference [29]. Two functions are applied to represent the uncertainty of sonar 
as follows: 


r(@) = 1—21(25)?, if 0 < |o| < 12.52, high 
0, if [6] > 12.52. 
P(p) = 1—(1 + tanh(2(p — py)))/2, (14.12) 


where @ represents the angle between the center-axis and the spot (i,j) measured in Fig. 
14.2. py is the pre-defined value, which reflects the smooth transferring point from the certainty 
to uncertainty. ['(@) shows that the nearer by center- axis is the spot (7,7), the larger is the 
density of the wave. I'(~) shows that the farther away from the sonar is it, the lower is the 
reliability, while the nearer by the sonar it is, the higher is the reliability of correct measurement. 


1) Probability Theory 


Elfes and Moravec [8, 9] firstly represented the probability of the grid occupied by obsta- 
cles with probability theory. Then Thrun, Fox and Burgard [27], Olson [17], Romero and 
Morales [19] also proposed the different methods of map reconstruction by themselves based on 
probability theory. According to the above methods, we give the general description of map 
building based on probability theory. To avoid an amount of computation, we suppose that 
all grids are independent. For every grid G;;, let s(Gi;) = E represent the grid empty, while 
s(Gj;) = O represent the grid occupied and P[s(Gi;) = E] and P[s(Gi;) = O] the probabilities 
of these events with the constraint P[s(Gj;) = E] + Pis(Gij) = O] = 1. According to the 
physical characteristics, the probability model to map the sonar perception datum is given by 
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(hey SPO aap ee 

[’1—-N (1-(2+4)")]/2, ifR-2<p<R-e, 
[1+ (1—a’)]/2, ifR-e<p<Rte, 
1/2, ifp>Rt+e 


P [s(Gij) = O|R] = P[s(o, 0) = O|R] = 


(14.13) 

where \’ = I'(0)-T'(p) and a = (p— R)/e. Seen from Eq. (14.13), the mapping between the 

sonar data and probability answers for the physical characteristics of sonar. For the data out- 

side the measurement range, the probability value is 0.5, that is, the uncertainty is the largest. 

When the distance between the grid in the range and the sonar is less than the measurement, 

the nearer by the sonar it is, the less is the possibility of grid occupied, while the nearer by the 
location of the measurement is the grid, and the larger is the possibility of grid occupied. 


The fusion algorithm for multi-sensors is given according to the Bayesian estimate as follows: 


P[s(Gij) = O|Re+i] - P [s(Gij) = O|Pi,- +> , Re 


~~ Pls(Giz) = X| Resi] - P[s(Giy) = X|R1,--- , Re] 
XE{E,O} 


P [s(Gij) = O|R1,--- , Resi] = 


(14.14) 
Remark: In order to make the equation (14.14) hold, we must suppose at the beginning of 
map building 
VGij EU, P [s(Gi;) = E| =P [s(Gi;) — O} = 0.5 


2) Fuzzy Set Theory 


Map building based on fuzzy logic is firstly proposed by Giuseppe Oriolo et al. [18]; they 
define two fuzzy sets V (represents grid empty) and 2 (represents grid occupied), which of size 
all are equal to U, correspondingly, their membership functions are wy and pq. Similarly, we 
can get fuzzy model to map the sonar perception datum. 


pS POY(Giz) =X - fu(p, R) (14.15) 


wPO(Gi;) =X'- fale, R), (14.16) 


where 


ke, if0<p<R-e, 
fu(p,R) = ke((R-p)/e)?, ifR-e<p<R, 
0, ifp>R. 


ko, if0<p<R-e, 
fale, R) = ko((R — p)/e)?, Lh-eapa Re; 
0, ifp>R. 
Here, fu(, R) represents the influence of p, R on the membership py of grid Gij. fo(p, R) 


represents the influence of p,R on the membership pig of grid Gj;. ke and ko are constants 
with 0 < kg < 1land0< ko <1. % is same as the definition in (14.13). Seen from the Eq. 
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(14.15) and (14.16), if the grid has the great possibility of occupation, then the membership is 
pow small, while the membership pug is great. If the grid is outside the measuring range, then 
bw = 0 and pe = 0. The fusion algorithm for multi-sensors is given according to the union 
operator in the fuzzy theory as follow: VX € {V,Q}, 


pee Gag) = wy (Gig) tg PO) — pe) Gy) pe” (Gig). (14.17) 


Remark: Initially, we suppose pry (Gi;) = ba(Gi;) = 0 (VGi; € U). According to membership 
of every grid, we can get the final map representation as follows: 


M=W2nN9QNANT (14.18) 
Here, A= UNO, 1 =VNQ. At the same time, the rule of fuzzy intersection is 
Hing (Gig) = ba(Giz) + Hy (Gig), V Gig € U, 
and the rule of fuzzy complementation is 
b(Giz) = 1 — pai(Giz), V Gig € U. 


The larger is the membership of Gi; belonging to the fuzzy set M, the greater is the possibility 
of the grid occupied. 


3) Dempster-Shafer Theory (DST) 


The idea of using belief functions for representing someone’s subjective feeling of uncertainty 
was first proposed by Shafer [21], following the seminal work of Dempster [4] about upper and 
lower probabilities induced by a multi-valued mappings. The use of belief functions as an alter- 
native to subjective probabilities for representing uncertainty was later justified axiomatically 
by Smets [26], who introduced the Transferable Belief Model (TBM), providing a clear and 
coherent interpretation of the various concepts underlying the theory. 

Let 6; (¢ = 1,2,...,n) be some exhaustive and exclusive elements (hypotheses) of interest 
taking on values in a finite discrete set 0, called the frame of discernment. Let us assume that 
an agent entertains beliefs concerning the value of 6;, given a certain evidential corpus. We 
postulate that these beliefs may be represented by a belief structure (or belief assignment), 
ie. a function from 2° to [0,1] verifying > 4¢@m(A) = 1 and m(0) = 0 for all A C 0, the 
quantity m(A) represents the mass of belief allocated to proposition ”6; C A”, and that cannot 
be allocated to any strict sub-proposition because of lack of evidence. The subsets A of © such 
that m(A) > 0 are called the focal elements of m. The information contained in the belief 
structure may be equivalently represented as a belief function bel, or as a plausibility function 
pl, defined respectively as bel(A) = d)pc4m(B) and pl(A) = di paazg M(B). The quantity 
bel(A), called the belief in A, is interpreted as the total degree of belief in A (ie. in the 
proposition ”6; C A”), whereas pl(A) denotes plausibility of A, i.e. the amount of belief that 
could potentially be transferred to A, taking into account the evidence that does not contradict 
that hypothesis. 

Now we assume the simplest situation that two distinct pieces of evidence induce two belief 
structures m, and mg. The orthogonal sum of m; and mg, denoted as m = m1 ® mz is defined 
as: 


m(A)=K7' S° mi(B)m2(C), (14.19) 
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i=D,let Be? (2) 


i-th spot of acquiring information 


Get sensors’ datum Robot’s location in global 
coordinate 


Is the grid in the intersection where 
sonar can scan? 


Is fusion times of this grid<=5 ? Stop fusing it 


I] y 
Go on fusing it, check the constraint 
and redistribute the conflict mass 


| check: Is the fusing the first times? Compute Bel 


Replace prior gbba of grid Is Bel of all grids 
with mass after fusion computed? 
N ag 


Exit 


Fig.9: Flowchart of procedure of sonar map building based on DSmT 


Figure 14.9: Flowchat of procedure of sonar map building based on DSmT and PCR5 


Here K = 1— SO mi(B)m2(C) for A 4 O and m(0) = 0. The orthogonal sum (also 


called Dempster’s rule of combination) is commutative and associative. It plays a fundamental 
operation for combining different evidential sources in evidence theory. DST as an information 
fusion method has been applied to the environment exploration and map reconstruction [1, 13]. 
This method can assure to have a precise result in fusing the same or different multi-sources 
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information from the sensors on the robot. According to the requirement of sonar grid map 
building, we need only to consider a simple 2D frame 0 = {61,42}, where 6; means ” the grid is 
empty”, and 9) means ”the grid is occupied”, and then we work with basic belief assignments 
defined on its power set 2° = {0,01, 02,0; U6}. According to DST, let mpgr(0) = 0, here 
mpst(91) is defined as the basic belief assignment function (bbaf) for grid-empty, mpgr(02) is 
defined as the basic belief assignment function (bbaf) for grid-occupied, mps7(41 U2) is defined 
as the basic belief assignment function for grid-ignorance. We may also construct basic belief 
assignment function such as mpsr(@1) = m(@1), Mps7r(42) = m(42), Mps7T(A1U82) = m(O,UG2). 
bbaf reflects still the characteristics of uncertainty for sonar grip map building in Fig. 14.2. 
Though here we define the same bbaf as DSmT, considering the definition of DST must be 
satisfied, we must renormalize them while acquiring sonar grid information [13, 14]. The new 
basic belief assignment after fusing two evidence sources from the sonar range finders can be 
obtained by the combination rule in Eq. (14.19). 


14.5 Simulation Experiment 


The experiment consists in simulating the autonomous navigation of a virtual Pioneer IT Robot 
carrying 16 simulated sonar detectors in a 5000mmx5000mm square array with an unknown 
obstacle/object. The map building with sonar sensors on the mobile robot is done from the 
simulator of SRIsim (shown in Fig.14.10) of ActivMedia company and our self-developing exper- 
imental or simulation platform together. (shown in Fig. 14.11). Here the platform developed 
with the tool software of visual c++ 6.0 and OpenGL servers as a client end, which can connect 
the server end (also developed by ourselves, which connects the SRIsim and the client). When 
the virtual robot runs in the virtual environment, the server end can collect many information 
(i.e. the location of robot, sensors reading, velocity .etc) from the SRIsim. Through the proto- 
col of TCP/IP, the client end can get any information from the server end and fuse them. The 
Pioneer II Robot may begin to run at arbitrary location; here we choose the location (1500mm, 
2700mm) with an 88 degrees angle the robot faces to. We let the robot move at speeds of 
transverse velocity 100mm/s and turning-velocity 50degree/s around the object in the world 
map plotted by the Mapper (a simple plotting software), which is opened in the SRIsim shown 
in Fig. 14.10. 


We adopt grid method to build map. Here we assume that all the sonar sensors have the 
same reliability. The global environment is divided into 50 x 50 lattices (which of size are same). 
The object in Fig. 14.10 is taken as a regular rectangular box. When the virtual robot runs 
around the object, through its sonar sensors, we can clearly recognize the object and know its 
appearance, and even its location in the environment. To describe the experiment clearly, the 
flowchart of procedure of sonar map is given in Fig. 14.9. The main steps of procedure based 
on the new tool are given as follows: 


1. Initialize the parameters of robot (i.e. initial location, moving velocity, etc.). 


2. Acquire 16 sonar readings, and robot’s location, when the robot is running (Here we set 
the first timer, of which interval is 100 ms.). 


3. Compute gbbaf of the fan-form area detected by each sonar sensor. 
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Figure 14.10: A world map opened in the SRIsim 


oF OR FERN S4F eth Oc ERM He FA TOA Seer eee 


nei send, 


wr] 
a) teen 
ern Mate vectes) pee | 


ae \eera 


os372 
—e 


Figure 14.11: The platform for simulation or real experiment 


4. Check whether some grids are scanned more than 5 times by sonar sensors (Same sonar in 
different location, or different sonar sensors. Of course, here we suppose each sonar sensor 
has the same characteristics.)? If ”yes”, then go to next step, otherwise, go to step 2. 


5. According to the combination rule and the PCR5 in (14.3) and (14.5) respectively, we can 
get the new basic belief masses, and redistribute the conflicting mass to the new basic 
belief masses in the order of the sequential fusion, until all 5 times are over. 
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6. Compute the credibility/total belief of occupancy bel(@2) of some grids, which have been 
fused according to (14.1). 


7. Update the map of the environment. (Here we set the second timer, of which interval is 
100 ms) Check whether all the grids have been fused? If” yes”, then stop robot and exit. 
Otherwise, go to step 2. 


Figure 14.12: Map reconstruction with DSmT coupling with PCR5 (3D) 


Finally, we rebuild the map shown in the Fig. 14.12 (3D) and Fig. 14.13 (2D) with the new 
tool. We also rebuild the map by other methods (i.e. probability theory, fuzzy theory, DST) in 
Fig. 14.14-Fig. 14.16. Because the other methods are not new, here we don’t give the detailed 
steps. If the reader has some interest in them, please refer to their corresponding reference. We 
give the result of comparison in Table 14.5. Through the figure and table, it can be concluded 
that: 

1) In Fig. 14.12, the Z axis shows the Belief of every grid occupied. The value 0 indicates that 
the grid is fully empty, and the value 1 indicates that this grid is fully occupied. This facilitates 
very much the development of human-computer interface of mobile robot exploring unknown, 
dangerous and sightless area. 

2) Low coupling. Even if there are many objects in grid map, but there occurs no phenomenon 
of the apparently servered, but actually connected. Thus, it supplies with a powerful evidence 
for self-localization, path planning and navigation of mobile robot. 

3) High validity of calculation. The fusion machine considering the restrained spreading arith- 
metic is adopted [29], and overcomes the shortcoming that the global grids in map must be 
reckoned once for sonar scanning every time, and improves the validity of calculation. 
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Figure 14.13: Map reconstruction with DSmT coupling with PCR5 (2D) 
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Figure 14.14: Map reconstruction with probability theory (2D) 


4) Seen from the Fig. 14.12, the new tool has a better performance than just DSmT in building 
map, (see also [13, 14]), because of considering the conflict factor, and redistributing the conflict 
masses to other basic belief masses according to the PCR5. 

5) Seen from the Table 14.5, Probability theory spends the least time, while Fuzzy theory spends 
the most time. But Map reconstruction with probability theory has very low precision and high 
mistaken judging rate shown in Fig. 14.14. Though the new tool spends a little more time than 
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Figure 14.15: Map reconstruction with fuzzy theory (2D) 
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Figure 14.16: Map reconstruction with DS theory (2D) 


probability theory. However, it has very high precision shown in Fig. 14.12 and Fig. 14.13 the 
same as the fuzzy theory and very low mistaken judging rate shown in Fig. 14.15. In fact, the 
comparison in map building between DST and DSmT have been made in details by us in [13]. 
Of course, here DST presents high precision and general mistaken rate shown in Fig. 14.16 
without considering the PCR rules and other conflict redistribution rules. We don’t compare 
the new tool from the fusion machine based on DST coupling with PCR5. Through the analysis 
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rating theory [Fuzzy Set Theory 
Total spent time (ms) 18146 35823 19673 19016 


Very high Very high 
Mistaken rate 


Table 14.1: Comparison of the performances of the different approaches 


of comparison among the four tools in Table 14.5, we testify the new tool to play a better role 
in map building. 


14.6 Conclusion 


In this chapter, we have applied a fusion machine based on DSmT coupled with PCR5 for mobile 
robot’s map building in a small environment. Then we have established the belief model for 
sonar grid map, and constructed the generalized basic belief assignment function. Through the 
simulation experiment, we also have compared the new tool with the other methods, and got 
much better performances for robot’s map building. Since it is necessary to consider also robot’s 
self-localization as soon as the size of environment becomes very large, complex, and irregular, 
we are currently doing some research works in Self-Localization And Mapping (SLAM) based 
on this new tool which improves the robustness and practicability of the fusion processing. In 
conclusion, our study has supplied a shortcut for human-computer interface for mobile robot 
exploring unknown environment and has established a firm foundation for the study of dynamic 
unknown environment and multi-robots’ building map and SLAM together. 
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